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I discuss the current theoretical interpretations of the anomalous J / ip suppression and 
their comparison to the most recent experimental data. 

1. Introduction 

The behavior of J/ip mesons in a hot strongly interacting medium was proposed as a 
test for its confinement status [Hj: it was argued that the J ft/), due to its small size and 
strong binding energy, can not break up as a consequence of interactions with normal 
hadrons, while in a deconfined medium, the color screening dissolves the cc bond. 

After this proposal, the study of the charmonium suppression in heavy ion collisions has 
aroused a lot of interest. With a careful and extensive analysis of the experimental data for 
different interacting systems (from proton-proton to proton-nucleus and nucleus-nucleus) 
it has been possible to observe a normal suppression of the J / ip meson, presumably due 
to the absorption of the preresonant cc state in the nuclear medium, in all interactions 
up to S-U collisions and peripheral Pb-Pb. The experimental data, for proton-nucleus 
(p-A) collisions, can be accurately described by a simple probabilistic formula for the J / ip 
survival probability 

S p a = ~^f- = J d 2 b J dzp A (b,z)exp^-(A-l) J dz' p A (b, z')a abs J , (1) 

where a abs is the effective cross-section for the absorption of the J/ip in nuclear matter: 
the experimental data fit gives a a b s = 4.18±0.35 mb [Ej; pA is the nuclear profile function, 
with parameters tabulated in [0j. 

The obvious generalization of Eq. Qto the nucleus-nucleus (A-B) case [E] gives an ex- 
cellent description of the J / ip suppression in S-U as well as in peripheral Pb-Pb collisions. 

In central Pb-Pb collisions, i.e. with impact parameter b < 8 — 8.5 fm, a stronger 
suppression is observed [ |3] : this anomalous suppression is the candidate as a signal of 
deconfined matter production at SPS. Recent NA60 data in In-In collisions confirmed the 
departure from the normal absorption in central collisions, with a pattern very similar to 
the one observed in Pb-Pb [ 5J. Surprisingly, also the first data from PHENIX (RHIC) 
show a suppression of the J/ip meson comparable with what observed at SPS energies[ 
U\ , while the common expectation pointed to a much stronger one. 
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The aim of this contribution is to give a critical overview of different theoretical inter- 
pretations of the anomalous suppression, comparing them to the most recent experimental 
data. In Sectional I will present a popular model which explains the NA38/NA50 data 
without assuming the deconfmement transition. Sections 0] and 03 present two different 
models of J/tp suppression in a deconfmed medium: thermal dissociation and parton per- 
colation. In Section El I will discuss the regeneration model: a recent approach which 
assumes that the observed J/ip's can be produced not only in the initial interactions but 
also in a subsequent stage of the system evolution. This second contribution is expected 
to be marginal at SPS, but very important at RHIC and LHC energies. Finally, in Sec- 
tion I will draw some conclusions of what we learned in these first ten years of J / ip 
suppression and what we need from future experiments. 

2. General remarks 

It is known from proton-nucleon and pion-nucleon [ |H] (and, more recently, from 
electron-proton interactions that a large fraction (about 30%-40%) of the observed 
J/ip's are the decay products of higher excited states of cc pairs (■*/>', x)- Since the life- 
time of these quarkonium states is much larger than the typical life-time of the medium 
which is produced in the nucleus-nucleus collisions, they decay in the vacuum. Therefore 
this medium (either hadronic gas or quark-gluon plasma) sees not only the ground state 
quarkonium, but also the different excited states, which have different properties (size, 
binding energies) and different behavior: a smaller binding energy (and, consequently, a 
larger radius) requires a lower dissolution temperature, in the case of a deconfmed medium; 
on the other hand, for the case of a hadronic system, a weakly bound quarkonium state 
has a large break-up cross-section for interactions with the other particles. 

Therefore the final J/ip survival probability, to be compared to the experimental data, 
has to be calculated as an average over the different components, each of them weighted 
with the corresponding fraction /j of contribution to the observed J/i/)'s in the final state: 

This fact has a very important consequence on the pattern of the J/ijj suppression as 
a function of the centrality and of the energy collisions and it must be considered for a 
careful comparison to experimental data. 

3. Suppression by hadronic interactions 

Before concluding that the anomalous suppression indicates the formation of deconfmed 
matter, one should check that the same affect can not be reproduced in a normal hadronic 
scenario. 

Many models of comover suppression have been proposed in the last years, the crucial 
parameter is the J/^-hadron inelastic cross-section a^h- The current theoretical estimates 
are quite vague, ranging from 0.1 to a few mb ! 

Since in S-U interactions the observed suppression is compatible with the normal one, 
the room for an additional suppression, due to inelastic interactions with secondary 
hadrons, is very small : this favors a very small a^h, since a large cross section would 
have sizable effects. 
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On the other hand, the anomalous suppression in central Pb-Pb is very strong and sets 
in quite abruptly : to obtain this effect a comover model must implement a mechanism 
which gives a very strong increase of the comover density in going from S-U to Pb-Pb 2 . 
The Dual-Parton Model (DPM) [E] is the most effective one in this respect. The num- 
ber of secondary hadrons is the sum of two contributions, proportional to the number of 
participating nucleons and to the number of individual nucleon-nucleon collisions respec- 
tively: 

N com (b, s, y, y/s) = Cx(y, y/s)N part (b, s, y/s) + C 2 (y, \fs)N coU (b, s, \fs). (3) 

In the previous formula, b is the impact parameter, s is the 2-dimensional coordinate in 
the transverse plane, y is the rapidity and i/i is the incident energy in the center of mass 
system. The coefficients C\ and C2 are predicted in the DPM, as well as their rapidity 
and energy dependence. The inelastic cross-section a^h is fitted to the experimental data 
and found to be 0.65 mb [H2]- Predictions for RHIC energies are, of course, possible [FTcA. 
with a correction taking into account shadowing effects (needed to reproduce the data on 
hadron multiplicity). 

In ref. [ ^] a reasonable description of S-U and Pb-Pb data is shown. However 
the model fails in describing the In-In results obtained by the NA60 Collab. [ 5j: the 
suppression predicted by the DPM is much stronger than what is observed, signaling that 
the effect of the collision term in EqElis too strong: one could try to reduce it by refitting 
the parameters, but then the agreement with Pb-Pb data would be lost. The same problem 
occurs at RHIC energies, where the experimental data are largely underestimated by the 
theoretical curve [EH- In my opinion, these data clearly disfavor the comover interactions 
as an explanation for the observed J/ip suppression. 

4. Thermal dissociation 

In sufficiently hot deconfined matter, color screening dissolves the binding of the quark- 
antiquark pair, and a stronger binding energy requires a higher temperature to be dis- 
solved. On a microscopic level, it was argued that only a hot medium provides sufficiently 
hard gluons to dissociate the quark- ant iquark bound state, and this again implies a dis- 
sociation hierarchy as function of the binding energy. Another possible mechanism is the 
decay into open charm mesons due to in-medium modification of mesonic masses [PH]. 

Implicit, in this approach, is the assumption that the medium probed by the quarkonium 
state is in thermal equilibrium. Lattice studies show that the transition between confined 
and deconfinement medium occurs at the critical temperature T c ~ 150 — 200 MeV. 

The large values of the charm (and bottom) quark mass allows potential theory to 
describe quite accurately the quarkonium spectroscopy [ [13]. 

The authors of ref. [H3] use the heavy-quark potential V(T) calculated on the lattice[ 
Ho] to obtain the quarkonium binding energy Mq(T) (i.e. the mass) and radius tq(T), 

2 One should, of course, make sure that the comover density needed to reproduce the data is not unreal- 
istically large: this problem, unfortunately, is usually neglected. The authors of ref. use a different 
approach: they set a physically meaningful limit for the density and the temperature of the hadron gas 
and conclude that it is not possible to describe the data of the anomalous suppression in the hadronic 
scenario 
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as a function of the temperature, by solving the Schrodinger equation. Below the critical 
temperature T c one finds that the masses Mq(T) of the cc bound states are a decreasing 
function of the temperature: when the condition Mq(T) < V^iT) for a given bound state 
is satisfied (V^o(T) is the asymptotic limit of the potential for r — > oo) the dissociation 
process QQ — > Qq + Qq is energetically favored and that bound state dissolves. 

Above the critical temperature T c one compares the binding radius tq(T) to the screen- 
ing radius l//i(T) (the distance at which the interquark interaction vanishes, according 
to the lattice results): again it is natural to assume that the bound state dissolves when 
tq(T) > \/ji{T). With these ingredients the authors of ref. [ [T3] find that x an d i>' 
dissolve already below the critical temperature, while the J/ip dissolves slightly above 
T c . As a consequence, a characteristic sequential suppression pattern is found, with three 
steps for the onset of the suppression of the if;', \ and directly produced J/ip respectively. 
One has to assume a model to translate the temperature into some experimental variable 
related to the centrality of the collision, and to take into account fluctuations in the num- 
ber of participants and collisions for a given impact parameter: the final result will be a 
smooth curve with two evident drops, one for the x an d the ip' suppression (if they occur 
for peripheral collisions they can hardly be separated) and the second one for the J/ip. 

5. Parton percolation 

Hadrons are made by partons. When two ore more hadrons overlap, their partons 
interact. In a normal hadron-hadron interaction, the overlapping phase is so short in time 
that the partons separate again into hadrons before reaching an equilibrium condition. 
On the other hand, in a nucleus-nucleus collision, the number of interacting hadrons is so 
high that their partons can interact several times, they therefore loose their "identity", 
so that they do not belong anymore to a particular hadron but form a big cluster of 
deconfined medium: the quark-gluon plasma (QGP). 

The percolation theory is a mathematical tool which studies how simple objects form 
clusters; there are applications of the percolation idea in many physical problems and 
the deconfinement transition in strongly interacting matter is one of them. For instance, 
in ref. [E] hadrons interact by exchanging color strings. When many hadrons interact 
simultaneously in a small space-time region, these strings overlap and, when their density 
reaches a critical value, they percolate. The model of hadron interaction based on color 
string exchange (below the percolation threshold) is able to reproduce many features of 
experimental data (see references in [EI)- This model therefore interpolate nicely from 
interactions in a normal hadronic medium and deconfined matter. A similar approach 
was followed in [fTS]: a model of string fusion and percolation is used to describe several 
experimental observables, including J/ip suppression. 

The work of references [01201121] is essentially focused on J/ip suppression by parton 
percolation, inspired by lattice results where it was shown that the deconfined transition 
in SU(2) Gauge Theory can be described by percolation of Polyakov Loops [122]. 

Parton percolation is an essential prerequisite for QGP formation. It should be noted 
that thermal equilibrium is not required and this is main difference with respect to the 
approach described in section HJ 

Cluster formation in percolation theory shows critical behavior : the cluster size diverges 
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at the percolation onset determined by the critical density n c of overlapping objects. In a 
finite system the percolation onset is defined as the point at which the grow of the cluster 
is more rapid [CHI- 

The above considerations apply to a nucleus-nucleus collision if one makes the following 
assumptions : the overlapping objects, forming clusters, are colored partons and the two- 
dimensional space where they are distributed is the transverse plane projection of the 
overlapping region of the two incident nuclei. A parton cluster represents a region in 
which color charges can move freely: if it extends over the entire space, one has, by 
definition, color deconfinement. 

In high energy nuclear collision there is one additional difficulty: the partons are emitted 
by the interacting nucleons, therefore their distribution is not uniform in the transverse 
plane, but rather reflects the original distribution of the initial nucleons, with a higher 
concentration in the center than near the surface of the nucleus. It is possible that in 
a given collision, only the most central region of the produced medium is dense and hot 
enough to allow for the deconfinement transition. In this situation a local definition of 
percolation, as the one used in ref. [[25J, is more appropriated since the J /if) meson 
is very small and therefore it is sensitive to the properties of a small spatial region. 
The percolation approach, therefore, has this evident advantage: it naturally allows the 
deconfined matter to be formed only in a limited part of the produced medium. 

Having specified the objects which form clusters, the partons, one has to describe their 
distribution in space to apply the percolation idea to nuclear collisions. Since these partons 
are emitted by the incident nucleons as a consequence of the strong interactions during 
the first stages of the nuclear collision, it is reasonable to assume that their number 
and spatial distribution is determined by the participating nucleons from which they 
originate, as done in refs [ I20J 1211 ■ Therefore the density of partons is given by the 
product of the participating nucleon density n s (b,A) (which depends on the nucleus A 
and on the impact parameter of the collision) and the number of partons per nucleon 
dN q (x,Q 2 )/dy (the parton distribution function, known from deep inelastic scattering 
experiments). The fraction x of the nucleon momentum carried by the parton is related 
to the incident energy yfs by x = {k T /yfs) (at midrapidity) ; k T is the average transverse 
momentum of the parton and it is inversely proportional to its transverse size. With these 
ingredients, the authors of ref. [ |^] find that the critical cluster density is reached, in 
Pb-Pb collisions at SPS, at b ~ 8 fm and the corresponding value of the average transverse 
momentum of the partons is Q c ~ 0.7 GeV. The scales of the charmonium states \ c and 
ip', given by the inverse of their radii calculated in potential theory, are about 0.6 and 
0.5 GeV respectively: they are therefore dissociated at the onset of percolation. On the 
other hand, directly produced J/ip's have smaller radii, therefore the average transverse 
momentum of the deconfined partons must be at least 0.9-1.0 GeV to resolve them, so 
only a denser medium, produced in more central collisions, can dissociate them. 

The J /if} survival probability in nuclear collisions is then obtained, from the above 
considerations, by assuming that about 40% of J/ip's (those coming from Xc and if)' 
decays) produced inside the percolating cluster disappear at the onset of percolation; 
those formed outside the cluster, i.e. near the surface, are not affected. The remaining 
60% of J/^'s (the directly produced ones) survive until a cluster of hard enough partons 
is produced (b ~ 3 — 4 fm). For a realistic comparison to the experimental data one has, 
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of course, to take into account impact parameter fluctuations (see, for instance, [HI). The 
result is in good agreement with the experimental data provided by NA50 Collab. 

The same model has be used to predict the survival probability as function of the 
centrality in In-In collisions at SPS. One finds, in this case, that the percolation onset is 
reached in semi-central collisions (N part ~ 140), but the threshold for the dissociation of 
directly produced J/ip's is never reached. One therefore expects, in this model, to observe 
a suppression pattern with only one step. The recently presented NA60 data show that 
this is indeed the case: the suppression pattern has one step only and the amount of the 
decrease is in remarkable agreement with the prediction of the percolation model, but the 
onset is wrong. It seems therefore that this model implements an excellent geometrical 
description of the deconfined phase (the amount of the suppression reflects the fraction 
of space occupied by the deconfined system) but needs to be improved. 

6. Suppression and Regeneration 

In the previous models it was assumed that the J j ip are produced exclusively upon the 
first, hard interactions of the colliding nucleons. Recently, the idea that the J/ip could 
also be created at the hadronization transition, by combining the cc pairs present in the 
plasma [[2Hj has been proposed. This phenomenon is possible if the J/ip states survive 
above the transition temperature, as suggested by lattice and potential calculations [ 

men. 

The authors of refs. [123 I2H1 121] consider both mechanisms : the primordial J/ip 
yield is subjected to nuclear absorption followed by dissociation in the QGP phase. A 
thermal contribution from statistical recombination of c and c quarks at the hadronization 
transition is added and also the J/ip break-up due to interactions in the hadron gas 
is included. The agreement with the experimental data of J/ip suppression in Pb-Pb 
collisions at SPS is excellent. The result also shows that the regeneration process gives 
a small contribution at SPS energies, while it is the dominant one at RHIC in central 
Au-Au collisions. 

7. Discussion and Conclusions 

The recently presented experimental data allows us to make a few considerations. First, 
they show that comover models can not describe simultaneously the S-U, Pb-Pb and In-In 
data at SPS energies. Even at RHIC energy they fail because they predict a suppression 
much stronger of what is observed. Actually, every model proposed so far which does not 
include the statistical regeneration fails in the comparison to PHENIX data. It would 
be tempting, therefore, to conclude that RHIC data show evidence of the regeneration 
mechanism at work. However, this conclusion is, in my opinion, premature. First of all, it 
is difficult to understand, even in the regeneration scenario, how can the J/ip suppression 
both in Au-Au and Cu-Cu at different energies (^Jsnn = 62 and 200 GeV) be equal: 
perhaps a different explanation should be invoked. 

A closer look to the newest NA50 and NA60 data shows that the J/tp suppression, 
and more precisely the ratio of the number of measured J/ip's to the expected ones 
(considering also the normal suppression) is always larger than 0.6. It was only in the 
previous data analysis that a few experimental points descended below this value for 
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very central collisions, suggesting the presence a 'second drop' due to the suppression of 
directly produced J/^'s. 

The newest SPS data show no evidence of direct J/ip suppression : the experimental 
data could be explained, with EqJ21 by assuming SjL = 1. This is confirmed by the 
agreement of the result of the percolation model with the NA60 data in central In-In 
collisions : in this model only ip' and x are suppressed, because the J/ip onset is out of 
reach in this system. The fact that the suppression in central In-In is about 0.8 instead 
of 0.6 (as it should be if one removes totally the ip' and % components) is due to finite 
size effects: the deconfinement sets in only in the central, hottest region and only the ip' 
and x formed inside are affected. 

This observation is corroborated by recent lattice calculations [ |23], where the dissoci- 
ation temperature for the J/ip in the QGP is found to be very high, about 1.6 — 2T C . 
Similar results have been found in recent potential model calculations, as shown in refs. [ 
I2H1 122| • The practical consequence of this fact is that the energy density required to break 
up the directly produced J ftp's can be as high as 30 GeV/fm 3 , definitely out of reach at 
SPS and even in central Au-Au collisions at RHIC ! If this is the case, we should be able 
to see the suppression of the direct J/ip component only at LHC 

In this way it is easy to understand why the suppression observed at RHIC is quantita- 
tively similar to the one seen at SPS: since the ip' and x components are already totally 
removed (except for finite size effects) at SPS, they can not be more suppressed at RHIC. 
J/ip's are unaffected both at SPS and at RHIC. Obviously, more precise data from RHIC 
experiments are necessary to confirm this scenario. 

The reason why all theoretical models not including regeneration effects underestimate 
the data is simply that all of them assume that the direct J/ip component can be sup- 
pressed: partially at SPS and more intensely at RHIC. 

It seems, therefore, that RHIC data can be explained without invoking the regenera- 
tion of J/ip's. Of course, one should check this idea, by looking, for instance, to other 
observables related to the J/ip, like flow, rapidity and px distribution. The regeneration 
model has already provided quantitative predictions in this respect [ ISTilEJT] . 

Finally, I would like to present a few comments about the deconfinement models in 
general. It is very common to assume that all charmonium states inside the QGP bubble 
are inexorably and instantaneously suppressed. This is, of course, an extreme hypothesis: 
even in a deconfined phase, some time is required to break up a J/^ (or ip' or x) meson, 
at least the time needed for the c and c to fly apart from each other. If the QGP lifetime 
is shorter of this time or if the cc is produced near the QGP surface with a transverse 
momentum enough to escape before being affected, then the corresponding charmonium 
state has a chance to survive even in the deconfined medium ! These effects are reasonably 
very important in small system (like In-In) or in peripheral collisions, right at the onset 
of the deconfinement. It was reasonable, so far, to neglect these corrections to keep under 
control the number of free parameters of the theory, but the level of accuracy reached 
now by the latest data forces us to seriously think about this necessary step forward in 
future developments. 
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